A Nonsingular Cosmology with a Scale-Invariant Spectrum of Cosmological 

Perturbations from Lee- Wick Theory 
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We study the cosmology of a Lee- Wick type scalar field theory. First, we consider homogeneous 
and isotropic background solutions and find that they are nonsingular, leading to cosmological 
bounces. Next, we analyze the spectrum of cosmological perturbations which result from this 
model. Unless either the potential of the Lee- Wick theory or the initial conditions are finely tuned, 
it is impossible to obtain background solutions which have a sufficiently long period of inflation 
after the bounce. More interestingly, however, we find that in the generic non-inflationary bouncing 
cosmology, perturbations created from quantum vacuum fluctuations in the contracting phase have 
the correct form to lead to a scale-invariant spectrum of metric inhomogeneities in the expanding 
phase. Since the background is non-singular, the evolution of the fluctuations is defined unam- 
biguously through the bounce. We also analyze the evolution of fluctuations which emerge from 
thermal initial conditions in the contracting phase. The spectrum of gravitational waves stemming 
for quantum vacuum fluctuations in the contracting phase is also scale-invariant, and the tensor to 
scalar ratio is not suppressed. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Recently, ideas originally due to Lee and Wick [2| were 
used to propose a "Lee- Wick Standard Model", a 
modification of the Standard Model of particle physics 
in which the Higgs mass is stabilized against quadrat- 
ically divergent radiative corrections and which in this 
sense is an alternative to supersymmetry for solving the 
hierarchy problem. The Lagrangian includes new higher 
derivative operators for each field. These operators can 
be eliminated by introducing a set of auxiliary fields, one 
for each field of the original model. The higher derivative 
terms have opposite sign for both the kinetic and mass 
terms, which indicates how the quadratic divergences in 
the Higgs mass can be cancelled. 

Fields with opposite sign of the kinetic term in the 
action have recently been invoked in cosmology to pro- 
vide models for dark energy. Fields with negative kinetic 
energy but positive potential energy are called "phan- 
tom field" [3| and were introduced to provide a possible 
mechanism for obtaining an equation of state of dark en- 
ergy with an equation of state parameter w < —1, where 
w = p/ p, p and p being pressure and energy density, 
respectively. In addition to the conceptual problems of 
having phantom fields (see e.g Q), phantom dark en- 
ergy models lead to future singularities. To avoid these 
problems, the "quintom model" Q was introduced. This 
model contains two scalar fields, one of them with a reg- 
ular sign kinetic term, the second with an opposite sign 
kinetic term. This model allows for a crossing of the 
"phantom divide" , i.e. a transition of the equation of 
state from w < —1 to w > —1. When applied to early 
universe cosmology, quintom models can lead to nonsin- 
gular cosmological backgrounds which correspond to a 



bouncing universe [6| [70|. 

The Higgs sector of the Lee- Wick Standard Model has 
similarities with the Lagrangian of a quintom model: the 
Higgs field has regular sign kinetic term but the auxiliary 
field has a negative sign kinetic term. Thus, it is logical 
to expect that the Lee- Wick model might give rise to 
a cosmological bounce and thus solve the cosmological 
singularity problem, in addition to solving the hierarchy 
problem. In this paper we show that this expectation is 
indeed realized. 

Given that the Lee- Wick model leads to a cosmological 
bounce, the cosmology of the very early universe may be 
very different from what is obtained by studying the cos- 
mology of the Standard Model. It is possible to introduce 
a potential for the scalar field in order to obtain a suffi- 
ciently long period of inflation after the bounce in order 
to solve the problems of Standard Big Bang Cosmology 
and to obtain a spectrum of nearly scale-invariant cosmo- 
logical fluctuations. However, this requires fine-tuning of 
the potential. On the other hand, given a bouncing cos- 
mology it is possible that the cosmological fluctuations 
originate in the contracting phase, as in the Pre-Big-Bang 
[9] or Ekpyrotic ^ scenarios. In this paper, we study 
the generation and evolution of fluctuations in our Lee- 
Wick type model. We consider both vacuum and thermal 
initial conditions for the fluctuations in the contracting 
phase and follow the perturbations through the bounce, 
a process which can be done unambiguously since the 
bounce is non-singular. 

We flnd that initial quantum vacuum fluctuations in 
the contracting phase have the right spectrum to de- 
velop into a scale-invariant spectrum in the expanding 
phase. What is responsible for this result is the fact that 
there is a coupling of the growing mode in the contract- 
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ing phase to the dominant (constant in time) mode in 
the expanding phase, and that this couphng scales with 
co-moving wave-number as k^. The Lee- Wick model thus 
leads to a concrete realization of the proposal of [llj 
(see also 0, H, and more recently to obtain a 
scale-invariant spectrum of fluctuations from a matter- 
dominated contracting phase (see also (14 ] for an analysis 
of gravitational wave evolution in this background). 

The outline of this paper is as follows. In the follow- 
ing section we introduce the Lee- Wick scalar field model 
which we will study in the rest of the paper. In Section 
3 we study the background solutions of this model, tak- 
ing initial conditions in the contracting phase. We show 
that, at least at the level of homogeneous and isotropic 
cosmology, it is easy to obtain a bouncing cosmology. In 
Section 4 we study how cosmological fluctuations set up 
in the initial contracting phase pass through the bounce. 
The evolution of the fluctuations is well behaved, Section 
5 contains the computation of the spectrum of gravita- 
tional waves, starting from quantum vacuum fluctuations 
in the contracting phase. We end with a discussion of our 
results. 



II. A LEE- WICK SCALAR FIELD MODEL 

We will take our starting Lagrangian for the scalar field 
(p to be 

C ^ ^d,^d^^-^{d^r-lm^^-Vi$), (1) 

where m is the mass of the scalar field and V is its in- 
teraction potential. The second term on the right-hand 
side is the higher derivative term, involving a new mass 
scale M. 

As discussed in j by introducing an auxiliary field 
and redefining the "normal" scalar field as 

^ = $ + 4>, (2) 
the Lagrangian takes the form 

(3) 

We thus see that M is the mass of the new scalar degree 
of freedom, the "Lee- Wick scalar" which comes from the 
extra degrees of freedom of the higher derivative theory. 
Note that both the kinetic term and the mass term of 
the Lee- Wick scalar have the opposite sign compared the 
signs for a regular scalar field. One may worry that the 
theory is unstable because of the wrong sign of the kinetic 
term of the Lee- Wick scalar llTI I. However, as was 

argued in , the perturbative expansion can be defined 
in a consistent way and the theory is unitary. Building 
on these workds, a recent study shows that Lee- Wick 
electrodynamics can be defined consistently as a ghost- 
free, unitary and Lorentz invariant theory [19|. 



By rotating the field basis, the mass term can be diag- 
onalized. However, the coupling between the two fields 
in the interaction term remains. To be specific, we con- 
sider a quartic interaction term. Thus, the Lagrangian 
we study is 

(4) 



III. BACKGROUND COSMOLOGY 

In this section we study the background cosmologi- 
cal equations which follow from coupling the matter La- 
grangian (HI) to Einstein gravity. For a homogeneous, 
isotropic and spatially flat universe the metric of space- 
time is 

ds^ = dt^ - a{tfdx^ , (5) 

where t is physical time, and x denote the co-moving 
spatial coordinates. The system of equations of motion 
consists of the Klein- Gordon equations 

4) + 3H4> + M^^ = ~X{(j)^4>f (6) 

for the two scalar fields and the Einstein expansion equa- 
tion 

= ^[^^'-^^' + ^™V-^Af^^^ + ^(0-0)^], 

(7) 

where H — a /a is the Hubble expansion rate and G is 
Newton's gravitational constant. An overdot denotes the 
derivative with respect to t. Combining these equations 
leads to the following expression for the change in the 
Hubble expansion rate 

H ^ -47rG((/)2 ~ . (8) 

from which we immediately see that it is possible for 
the background cosmology to cross the "phantom divide" 
H = 0. 

Let us take a first look to how it is possible to obtain 
a bouncing cosmology in our model. We assume that 
the universe starts in a contracting phase and that the 
contribution of (j) in the equations of motion dominates 
over that of the Lee- Wick scalar. This will typically be 
the case at low energy densities and curvatures. As the 
universe contracts and the energy density increases, the 
relative importance of compared to 4> will grow. From 
([7|) it follows that there will be a time when H = - this 
is a necessary condition for the bounce point. From ([5]) 
it follows that at the bounce point H > 0. Hence, we 
indeed have a transition from a contracting phase to an 
expanding phase, i.e. a cosmological bounce. 
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Let us now consider the above argument in a bit more 
detail. For the moment we will set the interaction La- 
grangian to zero, i.e. we will assume A = 0. We begin the 
evolution during the contracting phase when the energy 
density is sufficiently low so that we expect the contribu- 
tion of the Lee- Wick scalar to the total energy density to 
be small. For these initial conditions, both matter fields 
will be oscillating, and the equation of state will hence be 
that of a matter dominated universe. In fact, as follows 
from the Klein- Gordon equations ^ which in this case 
reduce to 

4> + 3H4) + m,'^(t) = 

4> + 3H^ + AP4> = (9) 

both scalar fields will be performing oscillations with am- 
plitudes A{t) and A{t) which are blue-shifting (i.e. in- 
creasing) at the same rate 

A{t) ^ A{t) ^ a{t)-^^\ (10) 

Eventually, the oscillations of the field 4> will freeze out. 
From studies of chaotic inflation ^20] it is well know that 
this happens when the amplitude A becomes of the order 
of the Planck mass rripi , more specifically when 

A = (12^)-i/2mpz . (11) 

After freeze-out, (j) will slowly roll up the potential and 
the equation of state will shift from w = to w ~ ~1 
(where w = p/p, p and p denoting pressure and energy 
density, respectively) leading to a deflationary phase dur- 
ing which the scale factor is decreasing almost exponen- 
tially. This phase is the time reversal of a period of slow- 
roll inflation. However, during this period the Lee- Wick 
field is still oscillating with rapidly increasing ampli- 
tude. Hence, its contribution to the energy density will 
rapidly catch up to that of (p. 

Let us give a rough estimate of the duration of the 
deflationary phase. It will depend crucially on the ini- 
tial ratio of the energy density of the Lee- Wick scalar (p 
to that of the regular scalar 0. Let us denote this ra- 
tio by JF. In the absence of coupling between the two 
scalar fields, i.e. for A = 0, the ratio will be unchanged 
during the period of matter domination when both fields 
are oscillating. However, once (j) enters the slow-rolling 
phase, the amplitude A will increase exponentially ac- 
cording to (jlOp while that of 4> will remain virtually un- 
changed. Thus, the condition on the duration At of the 
deflationary phase is 

\H\At ^ N ^ ^logm-') . (12) 

Thus, to obtain a deflationary phase with > 50 (which 
in the expanding phase will correspond to a period of in- 
flation of sufficient length to solve the cosmological prob- 
lems of the Standard Big Bang Model) required severe 
fine-tuning of the initial conditions. As we will discuss 



below, this problem may be even worse if coupling be- 
tween the two fields is allowed. 

Once the contribution of the Lee- Wick scalar to the 
energy density catches up to that of the original scalar 
field, the deflationary phase will end and a cosmological 
bounce will occur. Note that once H = 0, the Lee- Wick 
scalar is still oscillating whereas (j) is slowly rolling. Thus, 
H > and we indeed have a transition from contraction 
to expansion. This is a behaviour which is not possible 
for Einstein gravity coupled to matter satisfying the weak 
energy condition. However, due to the negative sign of 
the kinetic term in the Lagrangian, the weak energy con- 
dition is violated in our model. Note that in bouncing 
cosmologies obtained in higher derivative gravity models 
such as |21| , it is the higher derivative gravitational terms 
which, when interpreted as matter, lead to a violation of 
the weak energy condition. 

The duration of the bounce can be estimated as fol- 
lows: The maximal amplitude iJ„i of \H\ before and after 
the bounce is set by 

H,n ^ m, (13) 

since it is determined by the potential energy at the field 
value where the slow rolling of 4> begins. The amplitude 
of H at the bounce, denoted by H^, can in turn be esti- 
mated by 

•2 

Hb - 47rG0 - ATTGm^mli - , (14) 

where in the first step we have used the fact that the 
kinetic energy of (p is negligible at the bounce, and in 
the second step the fact that the bounce is determined 
by having the same absolute value of energy densities of 
(j) and (p, and that the field value of </> at the bounce is 
about rupi . The bounce time Atf, can now be determined 
via 

HbAh = 2H^ . (15) 

This gives 

Atb . (16) 

Note from the above that the value of i?,„ is set by the 
mass of (j), not the mass of (j). Similarly, the bounce time 
is determined by the mass of the original scalar field and 
not of its Lee- Wick partner. 

After the bounce, the amplitude of the oscillations of 
the Lee- Wick scalar exponentially decreases while is 
now slowly rolling down the potential. This is a phase of 
infiation which is time-symmetric to the phase of defla- 
tion before the bounce. As we have seen, without flne- 
tuning of the initial contribution of the Lee- Wick scalar 
to the energy density, the period of inflation will be too 
small for inflation to solve the various problems of stan- 
dard cosmology which inflation was invented to solve (2^ 
(see also [23| ) (such as the horizon and flatness problems) . 

Let us add some comments on the effects of allowing 
a coupling between the two scalar fields. We expect that 
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this will lead to a gradual flow of energy between the 
regular scalar and the Lee- Wick scalar such that at an 
energy density corresponding to the scale of the Lee- Wick 
scalar, the energy density in the Lee- Wick scalar will be- 
gin to dominate. Thus, allowing for A 7^ will lead to a 
shorter deflationary phase and may completely eliminate 
the period of deflation. Complete elimination of the de- 
flationary phase will occur if the energy density in cf) is 
larger than at cj) = (127r)~^/^mp;, where G = m~l^. 
This is the case if (making use of (fTTj) 

M < ((127^)-l/2„^TOp,)l/^ (17) 

The approximate analytical analysis summarized 
above is supported by exact numerical results. We have 



solved the coupled equations of motion for the scale fac- 
tor and the two scalar fields (j) and numerically. Figure 
1 presents the results in the base of a non-interacting 
model. We plot the time evolution of the scalar field 
(j) (denoted (pi in the figure), its Lee- Wick partner (f> 
(denoted by (^2) and the equation of state parameter 
w. As is evident, there is a non-singular cosmological 
bounce, there is no deflationary phase, but the equation 
of state parameter w crosses the phantom divide around 
the bounce point. Note that the scalar fleld (j) stops os- 
cillating near the bounce, whereas the Lee- Wick scalar 
continues to oscillate and therefore increases in magni- 
tude by a large factor during the latter stages of the con- 
tracting phase (which is why we have plotted the time 
evolution of on two different scales). 



Note that in the non-interacting model, the bounce is 
symmetric. In Figure 2 we present the corresponding flg- 
ure in the case of an interacting model with the value of 
A chosen to be A = 1.64 x 10~^^. In this case, the bounce 
is clearly asymmetric. As a second major difference com- 
pared to the simulation of Figure 1, the ratio of masses 
was chosen to be almost 100 in this case as opposed to 
only 2 in the first simulation. Because of the large ratio 



of the masses (and the corresponding initial conditions 
for which the energy in (f> greatly dominates over than in 
(/>), the background evolution enters a brief deflationary 
phase at the end of the contraction phase. However, due 
to the presence of interactions, the energy density in 
does not come to dominate again right after the bounce 
and hence the period of inflation which would be the time 
reversal of the phase of delation is absent. 



Finally, in Figure 3 we plot the number of e-foldings 
of the deflationary phase as a function of the ratio of 
to p^, in the model without interactions between the two 
scalar fields. As predicted by our analytical approxima- 
tions, the scaling of the period of deflation as a function 
of the ratio is roughly logarithmic. 



IV. COSMOLOGICAL FLUCTUATIONS 

A. General considerations 

It is useful to first consider the space-time sketch ^ 
of our non-singular bouncing cosmology. We choose the 
bounce time to correspond to t = Q. Long before the 
bounce, the equation of state is that of matter. During 
this period, the Hubble radius is decreasing linearly and 
H < 1. At a time denoted —tu (in analogy with the 
notation in inflationary cosmology) there is a transition 
to a period of deflation during which the Hubble radius 
is constant. However, as argued in the previous 
section, this period will be of short duration and ends 
at a time —ti when a brief bouncing phase covering the 
time interval [— < t < U begins. During this period 
H > 0. After the bouncing phase there is a short period 
of inflation lasting from ti to 1^, after which the universe 



enters a matter-dominated expansion phase with H < 0. 

In Figure [3] we also plot the evolution of the phys- 
ical length corresponding to a fixed co-moving scale. 
This scale is the wavelength of the fluctuation mode k 
(k standing for the co-moving wavenumber) which we 
want to follow. The wavelength begins in the matter- 
dominated phase of contraction on sub-Hubble scale, ex- 
its the Hubble radius during this phase at a time which 
we denote — i//(fc), and re-enters the Hubble radius dur- 
ing the matter-dominated phase at the time tnik). 

Note that if the energy density at the bounce point is 
given by the scale 77 of Grand Uniflcation (77 ~ lO^^GeV), 
then the physical wavelength of a perturbation mode cor- 
responding to the current Hubble radius is of the order of 
1mm, i.e. in the far infrared. In this sense, the evolution 
of fluctuations in this bouncing cosmology is free of the 
trans-Planckian problem [2J, which effects the evo- 
lution of fluctuations in all inflationary models in which 
the period of inflation lasts more than about 70 c-foldings 
(this number assumes that the scale of inflation is of the 
order of Grand Unification). 

Since our bounce is non-singular, the computation of 
the evolution of fluctuations is free of the matching con- 
dition ambiguities which affect the study of fluctuations 
in singular bouncing cosmologies such as the Ekpyrotic 
scenario (see [1^ [23, IH, [2^ [sO, HH for some early pa- 
pers on the problem of matching fluctuations through the 
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FIG. 1: Evolution of the background fields tj), tj) and of the background equation of state parameter m in a non-interacting model 
as a function of cosmic time (horizontal axis). The background fields are plotted in dimensionless units by normalizing by the 
mass Mrec ~ W~^mpi. Similarly, the time axis is displayed in units of M^J.. The mass parameters m and M were chosen to 

be m = bMrec and M = lOM^ec. The initial conditions were = 1.74 x lO^Mrec, (t>i = 1.44 x IQ* M'^^^, (ja = 8.98M,.ec, 0i = 

-14.08Mr^ec- 



bounce in Ekpyrotic cosmology). 

There is another important difference in the study of 
cosmological fluctuations between non-singular bouncing 
cosmologies and the inflationary scenario. It is usually 
argued that the exponential expansion of space during 
inflation red-shifts any pre-existing matter and the re- 
lated matter fluctuations, leaving behind a vacuum mat- 
ter state. Thus, perturbations in this setup are quantum 
vacuum fluctuations [tT]. On the other hand, in a bounc- 
ing cosmology the fluctuations are set up at low densities 
and temperatures in the contracting phase. There is no 
mechanism that red-shifts initial classical fluctuations. 
Thus, there is no reason to prefer vacuum over thermal 



initial perturbations. In the following, we will consider 
both choices. 



B. Equations for cosmological perturbations 

We begin by writing the metric including cosmological 
fluctuations in longitudinal gauge, assuming that there 
is no anisotropic stress (see [32] for a comprehensive dis- 
cussion of the theory of cosmological perturbations and 
[33| for a briefer survey) 

ds2 ^ q(^)2 ^ 2$)d?72 „ (1 _ 2$)dx2] , (18) 
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FIG. 2: Evolution of the background fields 0, and of the background equation of state parameter w in a non-interacting model 
as a function of cosmic time (horizontal axis). The background fields are plotted in dimensionless units by normalizing by the 
mass Mrec = 10~^mpi. Similarly, the time axis is displayed in units of The mass parameters m and M were chosen 

to be m = 1.4 x lO'^Mrec and M = WMrec- The initial conditions were 4>i = -3.57 x lO^Mrec, 4>i = 5.56 x 10^ M^^c, 4>t = 
2.98 X lO-eMrec, 4>i = -1.39 x 10"®M^^ec- 




FIG. 3: Plot of the duration of the deflationary phase as a 
function of the ratio of energy densities of </!> and (horizontal 
axis). The duration (vertical axis) is shown in terms of the 
e-folding number of deflation. 



FIG. 4: A sketch of the evolution of scales in a bouncing 
universe. The horizontal axis is co-moving spatial coordinate, 
the vertical axis is conformal time. Plotted are the Hubble 
radius \Ti\~°° and the wavelength A of a fluctuations with 
comoving wavenumber k. 
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where $(x, i) is the generahzed Newtonian gravitational 
potential which represents the metric fluctuations. It is 
convenient to write the equations in terms of conformal 
time 77 defined via dt = a{t)dr]. 

The Einstein equations linearly expanded in $ lead to 
the following equation of motion for the Fourier mode of 
<f> with co-moving wave-number k: 



(19) 



where the derivative with respect to conformal time is 
denoted by a prime, Ti. = a' /a, and dcj) is the fluctuation 
in (j). In deriving this equation, we have assumed that 
the background is dominated by the field (p. This will be 
the case except at the bounce. 

In inflationary cosmology, it has proven to be conve- 
nient to use the variable C, the curvature fluctuation in 
co-moving gauge, which in terms of $ is given by 



c = $ 



n 



(20) 



In any eternally expanding universe in which 1 + w 7^ 
0, the variable C is conserved on super-Hubble scales in 
the absence of entropy fluctuations (s^, [H, [11] , since - 
neglecting terms of the order - the equation of motion 
(fT9)) for <& is equivalent to 



{l + w)C ^ 0. 



(21) 



When considering the quantum theory of cosmological 
perturbations, it is important to identify the fluctuation 
variable which has canonical kinetic term. It is with re- 
spect to this variable, commonly denoted by v, that the 
canonical commutation relations must be imposed (see 
[stI . [sst for the quantum theory of cosmological pertur- 
bations). It turns out that the variable is simply related 
to C 



V = zC, 



(22) 



where the background variable z is the following com- 
bination of the background metric and the background 
matter field (f> (for simplicity we are assuming here only 
one matter field) 



IT 



(23) 



If the equation of state is constant in time, then z{r]) is 
proportional to a{r/). 

The equation of motion for v is 



\v = 0, 



(24) 



On sub-Hubble scales, it follows from (1241) that v is per- 
forming harmonic oscillations as a function of conformal 



time. On the other hand, on super-Hubble scales v is 
frozen in and v{r]) ~ z{ri). 

In terms of the variable z, the relationship between the 
metric fiuctuation <& and the canonical field v takes the 
form [H 



"(-)' 



n 



(25) 



The variable C has proven to be a convenient variable to 
use in inflationary cosmology. It was therefore taken for 
granted that it would also be a useful variable in bouncing 
cosmologies, and that it would remain conserved between 
when the mode k exits the Hubble radius during the pe- 
riod of contraction at the time — tff(fc) and the time tnik) 
of re-entry in the expanding phase. In the context of 
singular bouncing cosmologies, the Hwang- Vishniac (39j 
(Deruelle-Mukhanov [l^l) matching conditions for fluctu- 
ations across the singularity lead to the conclusion that 
C should be conserved. However, as pointed out in [3]| . 
the applicability of these matching conditions is ques- 
tionable since the matching conditions are not satisfled 
by the background. 

Non-singular bouncing cosmologies do not require ad- 
hoc matching conditions - the fluctuations can be fol- 
lowed through the bounce (as long as their amplitude 
remains sufficiently small such that linear perturbation 
theory does not break down) . As has recently been shown 
in several examples of non-singular bounces, the equation 
of motion for C develops singularities around the bounce 
point [U (43, [4^ , whereas the equation of motion for <I> 
remains well defined. One of the reasons for the singular- 
ities in the equation of motion for ( is that the co-moving 
gauge has a singularity at a cosmological bounce. Thus, 
in the following we will follow the evolution of the fiuc- 
tuations in terms of 

If the initial fluctuations in the contracting phase are 
due to thermal matter, then the initial values of <I> and 
its derivative follow from the perturbations in the energy 
density of matter. If, on the other hand, we assume vac- 
uum initial fiuctuations, then the initial inhomogeneities 
are given in terms of the canonical variable v and the 
initial values of and $ must be induced from v via the 
relation (1251). 



C. General solutions 

Let us briefiy review the general solution of the equa- 
tion of motion (|19p for $ on super-Hubble scales. We 
will keep the discussion quite general in this subsection 
and assume that the equation of state parameter is given 
by some w. In this case, a{t) scales as 



a{t) ~ tP 



with 



P = 



3{l + w) 



(26) 



(27) 
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From the definition of conformal time rj it then follows 
that 



r, ^ t 



i-p 



(28) 



The condition for the Hubble radius crossing time 
tnik) for a mode with co- moving wave- number k is 



Hence 



aitHik))k-' = H-\tH{k)) ^tnik). (29) 



m{k) ^ k-\ (30) 



As is well known, one of the two modes of $ on super- 
Hubble scales is constant, whereas the other is decaying 
in an expanding universe and growing in a contracting 
one. Specifically we have (see e.g. [ill) 



$(fc,r;) = D{k)+S{k)7j- 



-2i/ 



where 



2iy 



5 + 3w 
l + 3w ' 



(31) 



(32) 



and where D{k) and S{k) are independent of time and 
carry the information about the spectra of the two modes. 
In the following, we will determine the spectra of these 
two modes for various thermal and vacuum initial condi- 
tions. 



D. Thermal initial conditions 

Here we assume that the initial fluctuations are given 
by thermal matter perturbations. As was done in the 
case of string gas cosmology [3, [i^ (see [i^ for a recent 
comprehensive review), we follow the matter perturba- 
tions up to Hubble radius crossing and then convert to 
metric fluctuations by making use of the perturbed Ein- 
stein constraint equation (the time-time component of 
the perturbed Einstein equations) which reads 



(33) 



In the above, STq is the fluctuation in the energy den- 
sity, and V is the co- moving spatial gradient. At Hubble 
radius crossing all three terms on the left-hand side of 
the above equation are of the same order of magnitude. 
Hence, modulo a constant of the order 1, the Fourier 
space correlation function of $ becomes 

< |$(fc)|2 > = 16^2g2^-V < |5T0(fc)|2 ^ (34) 

where the pointed brackets indicate ensemble averaging. 

The energy density fluctuations are determined by 
thermodynamics. First, we express the momentum space 
energy density correlation function for co-moving wave- 
number k in terms of the r.m.s. position space mass 
fluctuation 



5M{R)^ ^ < |<5T°(fc)|^ 



> 



(35) 



where R — ak~^ is the physical radius of the region cor- 
responding to the wave-number k. The mass fluctuations 
are determined by the specific heat capacity Cy 



SM{R)^ = T^CviR), 



(36) 



where T is the temperature of the system. For a gas of 
point particles, the heat capacity is proportional to R^, 



Cv{R) = cyr^i?^ 



(37) 



where cy is a constant. Note that this result is in agree- 
ment with the intuition that on scales larger than the 
thermal correlation length ^ the heat capacity scales 
as a random walk. 

Inserting and into ([M]) we obtain the 

following power spectrum for $ 



< |$(fc)|2 IGTr^G^fc-^T^cy 



(38) 



We need to evaluate this expression at Hubble radius 
crossing tnik) since we will be using the corresponding 
value as the initial condition for the evolution of $ on 
super-Hubble scales: 

fc3 < |$(fc)|2 > {tnik)) = 16Tr^G^k~^T'^{tH{k))cv . 

(39) 

We now use (l39l) to infer the power spectra of the S and 
D modes in the case of thermal matter initial conditions. 
We are assuming that the initial value of $ at Hubble 
radius crossing gets distributed equally among the two 
modes The power spectrum Poik) of the constant mode 
D is the same as that of $ at Hubble radius crossing 



(40) 



where the second exponent comes from making use of 
T{tHik)) ~ a-\tH{k)) ^ tH{k)-P ^ fcT^ . (41) 

The power spectrum of S is the spectrum of $ at Hub- 
ble radius crossing modulated by the factor rinik)'^''- 



Ps{k) ~ fc-l+T^-'^^ 



(42) 



In the example we are interested in, fiuctuations leave 
the Hubble radius is the matter epoch and hence p — 2/3 
and u; = 0. Thus, from the above we see that the spectra 
of D and S scale as 



Poik) 
Psik) 



(43) 



The spectrum of the D mode is extremely blue. The blue 
tilt is due to the thermal suppression of the spectrum 
at large wavelengths. It is also easy to understand why 
the spectrum of the S mode is less blue than that of 
the D mode: the S mode grows on super-Hubble scales, 
and large wavelength modes experience the growth for 
a longer period of time. Our calculation shows that the 
difference in growth on super-Hubble scales dominates 
over the thermal suppression of long wavelength modes. 
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E. Vacuum initial conditions 

Vacuum initial conditions are given in terms of the 
canonically normalized variable Vk being in its quantum 
vacuum state [32i] 



(44) 



Inserting this into (P5)) and making use of the fact that on 
sub-Hubble scales the derivative of the oscillating factor 
dominates over the derivative of other terms leads to the 
following initial conditions in terms of <f>: 



(45) 



i.e. a spectrum which is proportional to fc~^/2. The same 
conclusion can be reached [271] by starting with vacuum 
fluctuations in the matter field (j) and inserting the result 
into the equations expressing $ and $ in terms of the 
matter field. Making use of ((23)) to eliminate z and of the 
background Friedmann equation to eliminate (f> in favor 
of H we find the following result for the power spectrum 
of 



P$(fc,ry) 



1 

2^ 



(46) 



TT nipi 



Making use of the definition of z from ((23)) it follows 
that the time dependent terms in ((45)) scale as H~^. 
Thus, 



(47) 



which allows us to evaluate the result at Hubble radius 
crossing 



$fc(t//(fc)) ^ k- 



(48) 



The above result ((48)) allows us to compute the spectra 
of both D and S modes of $ on super-Hubble scales, 
assuming - as we did in the previous subsection - that 
^k{tH{k)) sources both modes equally: 



<^D{k) 
<^s{k) 



kr 
k- 



In the case we are interested in p = 2/3, w 
2z^ = 5 we obtain 



*D(fc) 
$s(fc) 



(49) 
(50) 

and 

(51) 
(52) 



Note that, as pointed out in the S~ mode leads 
to a scale-invariant spectrum of fluctuations of C in the 
contracting phase. 



This compares to the results obtained for Ekpyrotic 
type contraction [27l | where p = 0, w = oo and 2i^ = 1 
and therefore 



$D(fc) 
$s(fc) 



(53) 
(54) 



which leads to a scale-invariant spectrum for the S mode 
which is growing in the phase of contraction. 

As we expect from the Hwang- Vishniac (Deruelle- 
Mukhanov) matching conditions, the S~ mode will cou- 
ple with a k^ suppression to the dominant mode in the 
expanding phase. If this is realized, we will obtain a 
scale-invariant spectrum of curvature fluctuations in the 
expanding phase. In the following subsection we will 
evolve the fluctuations through the non-singular bounce 
and infer the spectrum at late times. We indeed find a 
late-time scale-invariant spectrum. 



Evolution of the fluctuations through the 
bounce 



Let us step back and write down the equation of motion 
for <i> in a slightly modified form (which is equivalent to 
(Unj) except that we allow for a general speed of sound Cs 
which is equal to 1 in our scalar field model) 



where 



^'l + 2aH<^i+k^ci<i>k 



H 



0, 



2HH 



(55) 



(56) 



Contracting phase 
In the contracting phase the equation ((55)) takes the 
form 



1 + 2i^c 

■q - VB^ 



with 



$1. -I- fc2c2$,. 



3wc 



0, 



2(1 + 3wc 



(57) 



(58) 



where the subscript "c" indicates that we are discussing 
the contracting phase. The general analytical solution is 

$fc = iv - VB-)-"' { k--'^D^J,^[csk{f^ ~ fiB-)] 



+r=5_J_,Jc,fc(r/-77B_)]| , (59) 

where the coefficients £)_ and S*- can be determined by 
the initial condition of the gravitational potential as de- 
scribed in the two previous subsections for different sets 
of initial conditions. In the above, rjB^ is a fixed time 
that corresponds to when the singular bounce would oc- 
cur if the universe were to remain matter-dominated. 
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Note that, when the wavelength of the perturbation is 
larger than Hubble radius with k ^ |7i|, the asymptoti- 
cal form of <i>fe can be written as 



n = D- 

where we define 



2-=r(i + i/,) 



(60) 



(61) 



As discussed in previous subsections, the D- mode is 
constant and the mode is growing in a contracting 
universe. Note from the definition of we have to lead- 
ing order in fc 



5 -I- 2>Wc 
3(1 -Hwc) 



(62) 



Thus, to this order, in the contracting phase Cfc is de- 
termined by the constant mode of which is sub- 
dominant. As discussed in detail in [^t'I the S mode does 
effect C when /c^ corrections to the solutions are taken 
into account. 

If we match the asymptotic form for (f> to the 
initial power spectrum of $ (see (|46p) at Hubble radius 
crossing and assume that the initial power is equally dis- 
tributed into the two modes, we obtain 



1 
2"^2 



(63) 



where the subscript H stands for the time of Hubble ra- 
dius crossing. 

Bouncing phase 

As we have shown in the section on the background dy- 
namics, the contribution of the higher derivative terms 
in the Lee- Wick model becomes more and more impor- 
tant as the universe contracts and will lead to a non- 
singular bounce. Thus, the universe will exit from the 
phase of matter-dominated contraction at some time tB-^ 
and then the EoS of the universe will cross —1 and fall 
to negative infinity rapidly. Correspondingly, the Hub- 
ble parameter reaches zero and leads to a bounce of the 
universe at the time ^b- After the bounce, the Lee- Wick 
field will recover its normal state with the higher deriva- 
tive terms rapidly decreasing in importance. 

It is rather complicated to solve the perturbation equa- 
tion directly from Eq. ([55]) . In order to solve the equa- 
tion analytically, we need to make some approximations 
to simplify it. Our approximation consists of choosing a 
convenient modelling of the Hubble parameter near the 
bounce of the form 



at 



(64) 



with some positive constant a which has dimensions of 
and whose magnitude is set by the microphysics of 
the bounce, in our case by the mass M of the Lee- Wick 
scalar. The time of the bounce was chosen to be i = 0. 



In this case, we can obtain an analytical form for the 
comoving Hubble parameter in the bouncing phase: 



n = 



y 



Uv- Vb) 



1 - 

12 



VB) 



(65) 



where as denotes the value of the scale factor at the 
bounce point r]B- 

Since the above parametrization should be valid only 
in the neighborhood of the bounce point, the quadratic 
and higher order terms of \i] — riB \ can be neglected. Con- 
sequently, the perturbation equation takes the following 
form 



(66) 



The solution of this equation can be written as 

I 1 



EkHi[y^{r) - tjb)] + Fk iFi[ 
X exp[-y{ri - T]B'f] , 



2' 2 



,y{v-VB) 



(67) 



which is constructed from the l-th Hermite polynomial 
and a confluent hypergeometric function with 



3 



22/ 



(68) 



and two undetermined coefficients Ek, F^. These two 
functions are linearly independent, and their asymptot- 
ical behaviors are mainly determined by the parameter 
I. 

When Cgfc^ ^ y, i.e. the wave-number of the mode 
is larger than the mass scale of the bounce, then both 
functions are oscillating. This case was already studied 
in Ref. 

However, in the current paper we are interested in the 
opposite limit, the limit in which the wavelength is much 
larger than the inverse mass scale of the bounce, i.e. the 
limit when 



(69) 



As we have argued in the section of our paper on the 
background evolution, the bounce takes place very fast 
and thus the condition will be satisfied for all wave- 
lengths we are interested in. In this case, we can expand 
the solution of the perturbation equation in a power se- 
ries in terms of ^{rj — tib)- Then the solution is given 
by 



with 



Fk + Ek^iv - iib) + (-1 - l)Fkyiv - vb) 



Ek 
Fk 



21+' 



n-k) 



Ek 



r(¥ 



-Eh + Fk 



(70) 

(71) 
(72) 
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and the subscript "b" represents the bouncing phase. In 
this case, we have 



(73) 



Therefore, the conservation of Ck is reahzed by the mode 
Fk when the bounce is fast enough. 

Now we study how to estabhsh the coefficients Ek and 
Fk- We need to use the Hwang- Vishniac [s^ (Deruelle- 
Mukhanov [i^) matching condition to hnk the fluctu- 
ations in contracting phase with those in the bounc- 
ing phase at the momentum 773 Note that since we 
are matching two contracting universes across a non- 
singular surface, the background satisfies the matching 
conditions, unlike the situation in the Ekpyrotic scenario 
with a singular bounce, Thus, it is justified to apply the 
matching conditions [31 1 . 

The matching conditions say that both <i>fc and 



Ck = Ck 



3 T-P-W 



(74) 



are continuous on the matching surface of constant en- 
ergy density. Taking use of matching conditions in the 
solutions ([50)1 and ([70)l . we can obtain the following re- 
lations: 



EkVviVB- - Vb) = -(i + 2l)<i>% - Q\b- 



Fk 



2l)n + Ck\B- 



(75) 



These relations show that the constant and growing 
modes of gravitational potential get mixed during the 
bounce. However, if we consider large wavelengths com- 
pared to the duration of the bounce, the second relation 
shows that is indeed conserved across the bounce. 
Expanding phase 
After the bounce, the higher derivative terms of Lee- 
Wick field rapidly decay. Therefore, a phase of matter- 
dominated expansion starts at the time 7]b+- In the ab- 
sence of interactions between the two scalar fields, the 
background cosmology will be time-symmetric about the 
bounce point. During the period after riB+ the back- 
ground evolution is the time reverse of the contracting 
phase. In the case A ^ an asymmetric bounce is pos- 
sible. To render our analysis more general, we assume 
that the equation of state in the expanding phase is We 



which could be different from that in the contracting 
stage which is Wc- 

The equation of motion for the gravitational potential 
is similar as Eq. (|57p but with the indexes 



5 + 3We 
2(1 -I- 3We 



(76) 



and "B+^' instead of i^c and "_B— ". Then the solution 
on super-Hubble scales take the form 



with 



VB+ = VB+ 



(r/ - f]B+Y'''' 
2 1 

1 -f 3We + 



(77) 
(78) 



The D+ mode of the gravitational potential is constant 
in time, as is the D_ mode in contracting phase. How- 
ever, the role of the S mode is very different. In the ex- 
panding phase S+ is the sub-dominant decreasing mode, 
whereas in the contracting phase S- is the dominant ex- 
panding mode. Therefore, the dominant mode of the 
curvature perturbation in the period of expansion is . 
As we will shown in the following, it inherits contribu- 
tions from both Z)_ and S- since these modes mix during 
the bounce. 

To determine the coefficients of the two modes in the 
expanding phase, we need to apply the matching condi- 
tion again, this time at the surface 775+- A straightfor- 
ward calculation yields 



Eky/v{riB+ - 
Fk 



vb) 



= -{\ + '^mi-c'k\B^ 
2i)n + Q\B+ ■ 



(79) 



Note again that it is justified to apply the matching con- 
ditions since the universe is expanding on both sides of 
the matching surface and thus the background also sat- 
isfies the matching conditions. 

By combining Eqs. ([75]) and ([79]) , we can establish the 
relation between the gravitational potentials in contract- 
ing and expanding phases. Since 5+ is a decaying mode, 
we will not write down its expression and focus our at- 
tention instead on the dominant mode In terms of 
the modes in the contracting phase, it is given by 



^ ^ (5 + 3uO(l±^^_ 
+ (l + M;,)(5 + 3we) 



J 



3(H-We)^2^2 X f 775+ -?7i3 - ^ - 2LI_ 



3(l + i«c) 



-M+( 



■nB--VB 3(1 -I- Wc) {r]B--fiB-) 



M+D- +M-{D- 



7 +0(fc^) ' 



(80) 



where we defined the parameters 

2 1 



which are independent of k. 
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From the above result we see that both the con- 
stant and growing modes of gravitational potential in 
the contracting phase effect the dominant mode after 
the bounce. However, the growing mode is suppressed 
by fc^ on large scales whereas the constant one transfers 
through the bounce without a change in the spectral in- 
dex. These results agree with what is obtained using the 
matching conditions at a singular hypersurface between 
the contracting and the expanding phase, as shown in 

Inspecting our result (|80l). we see that there are two 
ways to obtain a scale- invariant spectrum of cosmological 
perturbations after the bounce. The first is to consider 
a model in which the D— mode in the contracting phase 
has a scale invariant spectrum, i.e. D ~ {k) ^ fc~^/^, the 
other is to take a scenario where S — (k) ^ fc~^/^. As 
follows from ([^ the first possibility is realized if p = oo, 
i.e. in an inflationary contracting phase. The second way 
is realized in the case of a matter-dominated contraction, 
a possibility already pointed out in fTl| (see also p7j). 
The Lee- Wick model yields a natural realization of this 
way. 

Let us now come back to our Lee- Wick background, 
and assume quantum vacuum fluctuations. We insert 
the values Wc = = 0, Cs = I into (|80p and assume a 
symmetric fast bounce. Thus 



5_ + 



5 24 



2fc2 



(82) 



As discussed in the subsection on vacuum initial condi- 
tions, then if the initial conditions are imposed at a time 
r] sufficiently early compared to the transition point -B— , 
we have oc fc~2 and therefore obtain Z)_ oc fca and 
S- oc k~'i . Substituting these relations into Eq. ([5^ . 
one can see that whereas the contribution of l)_ to the 
final spectrum oi D+ vanishes on large scales, the con- 
tribution of S- which starts out deep red is blue-tilted 
by exactly the right amount to yield a final spectrum 
proportional to k~2 which is the scale-invariant form. 

To compute the amplitude of the spectrum, we insert 
the values (|8ip and the expression (|63p for the value of 
S- into (j82p and use the background Friedmann equation 
to replace the Hubble parameter by the energy density. 
This yields 



10\/2Mp2 



(83) 



Therefore, the power spectrum of the gravitational po- 
tential in this case can be expressed as 

Note that as long as M <C Mpi, the power spectrum 
of metric fiuctuations remains much smaller than 1 and 
thus linear perturbation theory is applicable throughout 
the bouncing phase. 



G. Numerical analysis 

Our analytical calculations involve approximations. 
Specifically, in the contracting phase the scalar fields and 
hence the equation of state are oscillating. But in our 
analytical analysis we have replaced the time-dependent 
equation of state parameter by its temporal average. It 
is thus important to confirm the results by numerical in- 
tegration of the full equations, namely Eq. (|19p coupled 
to the equation for the scalar matter field fluctuation. 

At first sight, it appears that the equation con- 
tains a singularity at all turnaround points of (p. Such 
singularities are known from the study of the evolution 
of $ during reheating taking into account the oscillatory 
nature of the inflaton field [41, |4^. However, this singu- 
larity is actually not present. Let us consider in addition 
to the dynamical perturbed Einstein equation (|19p the 
perturbed Einstein constraint equation 

+ = 47rG(0'(5(/) - 4>'S4>) . (85) 

Inserting ((SSl) into (HH) yields 

$" + 67^$' + 2(7i' + 27^2)$-|-fc2$ 

= 87tG(2H+^)(I)'6(P, (86) 
(p 

from which it is clear that the singularity has disap- 
peared. Thus, we numerically solve (jBS]) coupled to the 
perturbed (j) equation 



2nS<j)'- 



,)S(t) = 4(/)'$' - 2a , (87) 



where the subscripts on V indicate the variables with 
respect to which the potential is differentiated. 

Figures and [5] show the results of our numerical in- 
tegration. The first figure shows the evolution in time of 
the metric fluctuation $ as a function of physical time 
(left side) and conformal time (right side) for different 
values of the comoving wavenumber k. We have chosen 
the bounce point to correspond to physical and confor- 
mal time 0. The initial conditions for $ were set at the 
initial time of the simulation according to the vacuum 
initial condition prescription discussed earlier. We see 
from this figure that before the bounce the perturbations 
are dominated by the growing mode. When the universe 
enters the bouncing phase, we see that the amplitude 
approaches a constant and passes smoothly through the 
bounce. The numerical evolution agrees well with the 
analytical solution (|70|) in the bouncing phase. After the 
bounce, the perturbations are dominated by the constant 
mode. The numerical evolution demonstrates that this 
mode can be inherited from the growing mode in con- 
tracting phase. 

Figure [5] shows the power spectrum of $ (lower panel) 
and the spectral index Us (upper panel) as a function of 
comoving wavenumber k. On large scales (small values of 
k), the power spectrum tends to a constant. The rise of 
the spectrum for large values of k is on scales which are 
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comparable to maximal value of the Hubble rate, i.e. for 
modes which have not spent time outside of the Hubble 
radius. 



10^ 
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FIG. 5: Result of the numerical evolution of the curvature 
perturbations with different comoving wavenumbers k in the 
Lee-Wick bounce. The horizontal axis in the left panel is 
cosmic time, and in the right panel it is comoving time. The 
initial values of the background parameters are the same as in 
Figure 1. The units of the time axis are M~2:, the comoving 
wavenumber k is unity for k = Mrec, as in Figure [2] 
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FIG. 6: Plot of the power spectrum of the curvature per- 
turbation $ (lower panel) and of the spectral index (upper 
panel) as functions of comoving wavenumbers k in the Lee- 
Wick bounce. The initial values of the background parameters 
are the same as in Figure 1. 



V. GRAVITATIONAL WAVES 

Now we turn to consider the evolution of gravitational 
waves (tensor perturbations) in our background, assum- 
ing they start out in the vacuum state on sub-Hubble 
scales in the contracting phase. Since at the level of lin- 
ear perturbation theory scalar metric fluctuations and 
gravitational waves decouple, we can focus on a metric 
containing only gravitational waves propagating in the 
background. The standard form of this metric in a spa- 
tially flat FRW background is 

ds^ = a{r]f[-dri^ + (% -1- h,^)dx'dx^] , (88) 

where the Latin indexes run over the spatial coordinates, 
and the tensor perturbation hij is real, transverse and 
traceless, i.e. 

hij = hji ; ha = ; /lyj = . (89) 

Due to these constraints, we only have two degrees of 
freedom in hij which correspond to two polarizations of 
gravitational waves. For each polarization state (labelled 
by r in the following), we can write hij(r],x) as a scalar 
field h^{ri, x) multiplied by a polarization tensor e^^ which 
is constant in space and time. 

If matter contains an anisotropic stress tensor dy , 
there is a non-vanishing source term in the equation of 
motion for tensor perturbations, namely 

n'- 

h'l^ + 2—h\j - V^/iy = 167rGaVy . (90) 

If matter consists of a set of canonically normalized scalar 
fields or a set of perfect fluids, there is no anisotropic 
stress and thus no source term at linear order in pertur- 
bation theory for gravitational waves. 

As usual, we go to Fourier space The Fourier trans- 
formations of the tensor perturbations and anisotropic 
stress tensor are give by, 

2 » ,3, 

%(r?,x) = ^ / h^{^^\,)el/^- . (91) 

In order to canonically quantize the gravitational waves, 
it is important to identify the variable in terms of which 
the action has canonical kinetic term. This variable turns 
out to be (see (HO] for a derivation) 

(where /i^ is a short hand notation for ft.*" (77, k)) in terms 
of which the Einstein action expanded to second order in 
becomes 

^ = il\ fim\'-{k'-^MWk. (93) 

r— 1 
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The resulting equation of motion for v'' is 

K)" + - —y = . 



(94) 



We are interested in computing the power spectrum 
of the tensor modes. Making use of ([5^ . it is related 
to the power spectrum of v (which is the same for each 
polarization state) via 



^AttG 



27r2 



The tensor spectral index ut is defined by 

dlnP^ 



riT 



dink 



(95) 



(96) 



The evolution of tensor perturbations is very similar 
to that of scalar perturbations. Initially the perturba- 
tions are inside the Hubble radius in the far past. Since 
the Hubble radius shrinks in the contracting phase, the 
modes with small comoving wave number exit the Hubble 
radius. After that the universe bounces to an expanding 
phase, so these Fourier modes will return into the Hubble 
radius. 

In the current paper we focus on a mode with small k 
so that it exits the Hubble radius in the contracting phase 
(rather than the bounce phase) , then passes through the 
bounce and finally re-enters the Hubble radius during the 
expanding phase. 

We divide the time interval into three periods like we 
did for the analysis of scalar metric fluctuations. Dur- 
ing the phase when the universe is contracting with an 
equation of state oscillating around w = 0, we have 



3[k{ri-fiB-)] 



VB-) 



(97) 



where f/B- — r/B- — 2/7Yb_. Here, the parameters 
and can be determined by the initial condition for 
gravitational waves, which is taken as the Bunch-Davies 
vacuum 



2k. 



(98) 



So we have 



Bl = - 



and 



(99) 



Therefore, the asymptotic form of the solution to the 
tensor perturbation in the contracting phase is 



During the bouncing phase, 
relation 



we have the approximate 



4 , 
—aa 

TT 



y_ 

3 



(101) 



Solving Eq. (f94|) . we have 

v{k,ri) = 

' cos [Z (77 



CI' 



3 ' 



Dl 



< 



(102) 



where we define l"^ 



|fc2 — ||. Since the Hubble pa- 
rameter approaches zero when the universe is transiting 
from the contracting to the expanding phase, all fluctu- 
ation modes return to the sub-Hubble region, but only 
for a very brief time. However, from the above solution 
we interestingly find that k^f^i^ k'^ /a\) and a) are 
comparable. 

After the bounce, an expanding phase with its EoS 
w = takes place. So the solution to the gravitational 
waves is given by 



= {-n-fiB+Y 



ElJ_s[k{r^-VB+)] 



+F^J,[k{rj^fjB+)] 



(103) 



where fB+ = ^54 
asymptotic form. 



-2/ni 



This solution takes on the 



2 3 



(104) 



V2 



k 2(ri-f)B-y 



outside Hubble radius; 



when the modes are outside the Hubble radius. 

Having obtained the solutions of the tensor pertur- 
bations in the different phases, now we need to match 
these solutions and determine the coefficients Cj, Z?|", 
E"! and respectively. This procedure is analogous to 
the matching process of scalar perturbations performed 
in the previous section. For a non-singular bounce sce- 
nario such as the bounce we are considering, the continu- 
ity of the background evolution implies that both v and 
v' are able to pass through the bounce smoothly. So we 
match V and v' in (|100p and (|102p on the surface r^-, 
and those in (|102p and (|103p on the surface tb+- With 
these matching conditions, we can determine all the co- 
efficients and finally get the solution for v at late times. 

Since in the specific model we considered in this paper, 
the evolution of the universe is symmetric with respective 
to the bounce point, we can simply take Hb- — —'Hb+- 
In addition, we have shown that the bounce takes place 
very fast on the time scale set by k~^ , so we have 
1{tb+ — tb^) S> 1. Therefore, we eventually obtain the 
approximate result 



— i=e ^'^^^ f)B-) inside Hubble radius. 



(100) 



/2k 



Fi 



'' 8A;3 



B + 



(105) 
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and the asymptotical form of v in the final stage can be 
expressed as 



(106) 



Now we are able to derive the power spectrum of pri- 
mordial gravitational waves. From the definition of Eq. 
(j95p . the primordial power spectrum is given by 



Prik) 



TT a 

'2'PB + 



277r2M4 



(107) 



From Eq. (|107p . we can read that the spectrum are scale- 
invariant on large scales (which is consistent with the 
result in Ref. 51]). 

Comparing our result of the tensor power spectrum 
with the result (|84p for the power spectrum of scalar met- 
ric fluctuations, we obtain a tensor to scalar ratio of the 
order of 30, which is in excess of the current observational 
bounds. The exact value of the ratio, however, will de- 
pend on the detailed modelling of the bounce phase (72| 
However, the conclusion that the tensor to scalar ratio 
will be rather large will be robust, and also agrees with 
the analysis of [47| done in a different context. 



VI. CONCLUSIONS AND DISCUSSION 

Recently, the Lee- Wick Standard Model has been sug- 
gested as an extension of the Standard Model of parti- 
cle physics providing an alternative to supersymmetry in 
terms of addressing the hierarchy problems. 

In this paper, we have considered the cosmology of the 
Higgs sector of the Lee- Wick Standard Model, an alterna- 
tive to supersymmetry to solving the hierarchy problem. 
We have found that homogeneous and isotropic solutions 
are non-singular. Thus, the Lee- Wick model provides a 
possible solution of the cosmological singularity problem. 

We then considered the spectrum of cosmological per- 
turbations and flnd that quantum vacuum fluctuations 
established in the contracting phase evolve into a scale- 
invariant spectrum in the expanding phase. Note that 
these results emerge without having to introduce any 
additional features into the model, unlike the situation 
in inflationary cosmology where the existence of a new 
scalar field satisfying slow-roll conditions must be as- 
sumed, or the situation in Ekpyrotic models where once 
again a scalar field with special features must be assumed. 

Tuning the amplitude of the spectrum of scalar metric 
fluctuations to agree with the amplitude inferred from 
CMB observations [s^], we can determine the scale AI of 
the new physics which is present in the Lee- Wick model. 
The required value of M turns out to be about lO^^GeV. 

We have also computed the spectrum of gravitational 
waves and also find a scale-invariant spectrum assuming 
that the fiuctuations are quantum vacuum in nature. The 



tensor to scalar ratio may be in excess of the current 
observational bounces, but the exact value will depend 
on the detailed modelling of the bounce phase. 

One of the main successes of cosmological inflation is 
the solution of the horizon, homogeneity, size and flatness 
problems of Standard Big Bang cosmology which it pro- 
vides. How does a bouncing cosmology such as our Lee- 
Wick model measure up against these successes? First 
of all, if the universe starts out large and cold, there are 
no horizon and size problems. If the spatial curvature at 
temperatures in the contracting phase comparable to the 
current temperature is not larger than the current spatial 
curvature, then there will be no flatness problem either 
because the deviation of Qk from decreases in the con- 
tracting phase at the same rate that it increases in the 
expanding phase. The key challenge for any bouncing 
cosmology is to control the magnitude of the inhomo- 
geneities and to provide a mechanism for preventing the 
universe to collapse into a gas of black holes at the end 
of the phase of contraction. For an attempt to address 
this issue in the case of string gas cosmology see (ssj . 

We would like to conclude this paper by putting our 
work in the context of previous work on perturbations in 
bouncing cosmologies. The issue of the mixing of the S 
and D modes at a cosmological bounce has been hotly 
debated in the literature since the Ekpyrotic scenario 
was proposed. In the case of the Ekpyrotic scenario, 
for vacuum initial conditions the S mode of inherits 
a scale-invariant spectrum, whereas the D mode obtains 
a blue spectrum with index n = 3 [1^, [27, 28, 29] (see 
also the more recent analysis of [HJ] and the recent re- 
view of [El]). This is also the spectrum of Accord- 
ing to the Hwang- Vishniac [s^ (Deruelle-Mukhanov [ioj ) 
matching conditions applied at a hypersurface on which 
we glue the expanding to the contracting universe, the 
mixing between the S— mode and the D+ mode is sup- 
pressed by a power of fc^ (see e.g. [28*] for a discussion 
of this point). Hence, the spectrum of metric fluctua- 
tions after the bounce is not scale-invariant. The Pre- 
Big-Bang scenario faces a similar problem [s^. These 
conclusions were confirmed in some specific models in 
which the bounce was smoothed out by making use of 
higher derivative gravity terms (see [53] in the case of 
the Pre-Big-Bang model and [H, [5^ in the case of the 
Ekpyrotic scenario). However, the use of the matching 
conditions was challenged in [31| where it was pointed 
out that if the background solution does not satisfy the 
matching conditions at the bounce, there is no reason 
to expect the fluctuations to do so. In fact, in the case 
of the Ekpyrotic scenario (which is intrinsically a higher- 
dimensional cosmology) , computations done in the higher 
dimensional framework yielded a successful transfer of 
the scale-invariant spectrum of metric fluctuations from 
the contracting to the expanding phase [13] , a conclusion 
which was confirmed in [6lj and, in a slightly different 
setting, in [12] [zl]. 

Calculations have also been done in some other non- 
singular bouncing models ^] . For example, studies done 
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in models in which the bounce is induced by a negative 
energy density scalar field found no unsuppressed match- 
ing between the growing perturbation mode in the con- 
tracting phase and the constant mode in the expanding 
phase [64l. l65|. in contrast to what was obtained in some 
initial work[3, Q. Both studies in models in which a 
bounce was generated by a curvature term in the Einstein 
action (66j and analyses in some other bouncing models 
[67l [68| yielded un-suppressed matching of the dominant 
modes of the contracting and expanding phases. 

The upshot of these analyses is that the transfer of 
fluctuations through a cosmological bounce can depend 
quite sensitively on the physics of the bounce. 

It was realized that the equation of motion for ^ has 
singularities in the case of a non-singular bounce, thus 
casting doubt on the belief that in all cases C is con- 
served at a bounce. It was shown that the $ equation 
is free of such singularities and is thus a safer equation 
to use [41I, im . In our previous work [i^ it was shown 
in the case of the quintom bounce model that there is 
unsuppressed mixing between the D+ and S— modes on 
length scales which are small compared to the duration 
of the bouncing period, whereas on longer length scales 
the mixing is suppressed (but not completely absent). In 
the present work, the bounce is short compared to the 
length scales we are interested in. 

Our work shows that the evolution of fluctuations 
through the non-singular bounce in the Lee- Wick model 
is rather standard. There is no un-suppressed coupling 
between the dominant modes of the contracting and ex- 
panding phases, and ^ is conserved at the bounce. 

In the current paper we have not considered radiation. 
Since the energy density in radiation increases faster than 
that in matter, radiation would dominate at early times. 
However, in the Lee- Wick standard model there is a Lee- 
Wick partner to each field. In particular, there is a Lee- 
Wick photon partner 7 of the radiation field 7. At high 
energy densities, then as a consequence of interactions 



between 7 and 7, we expect that energy will flow from 7 
into 7, like it flows from to (/> in our scalar field model. 
Then, a cosmological bounce would occur in a manner 
similar to how it occurs in our model. Adding inter- 
mediate phases of radiation between the bouncing phase 
and the contracting and expanding matter phases will 
not change our results concerning the spectrum of fluc- 
tuations for modes which exit the Hubble radius during 
the phase of matter domination, which are the modes we 
are interested in when trying to explain the large-scale 
structure of the universe and the CMB anisotropies. A 
study of these issues is left to a followup paper. 

It would also be interesting to consider entropy fluctu- 
ations and non-Gaussian signatures of our scenario. We 
leave these topics for future research. 

Note added: while this paper was being prepared for 
submission, a preprint appeared [6^ pointing out that 
the Lee- Wick model provides a realization of the quin- 
tom scenario and could be applied to study the current 
acceleration of the universe. We find it more natural 
to consider the corrections to the cosmological evolution 
which are obtained in the very early universe. 
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The use of a second scalar field to obtain a nonsingular 

bounce was already discussed in 0, @1 

In light of the trans-Planckian problem for inflationary 

fluctuations, one may view this prescription with some 

scepticism. 

We have modelled the contracting and expanding phases 
with a background with a constant equation of state. In 
this case, the squeezing factor for scalar and tensor cos- 
mological perturbations is the same. If the equation of 
state changes, then the amplitude of the scalar spectrum 
can be enhanced, as happens in inflationary cosmology at 
the time of reheating [35l.l36|. Since l-|-ui does not change 
by orders of magnitude in the Lee- Wick bouncing cos- 
mology - whereas it does change by orders of magnitude 
in the case of inflationary cosmology - we do not expect 
a dramatic enhancement of the scalar power spectrum. 
However, we do expect some increase in the amplitude. 
A detailed study of this issue is left to further work. 
It wold go beyond the scope of this discussion to review 
attempts to preserve a scale-invariant spectrum in Pre- 
Big-Bang and Ekpyrotic scenarios via introducing en- 
tropy fluctuations. 



